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Resonant e'^e^ pair production by an electron in a magnetic field near the process threshold has 
been analytically studied. Using the Nikishov's theorem an estimation of the number of events has 
been made in the magnetic field equivalent to laser wave in the SLAC experiment [D. Burke et al., 
Phys. Rev. Lett. 79, 1626 (1997)]. The obtained estimation is in reasonable agreement with the 
experimental data. 
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I. INTRODUCTION 

Fundamental processes in intense external electromag- 
netic fields are of great interest due to tlie existence of 
strongly magnetized neutron stars and the construction 
of high-power laser systems. Known physical processes 
are modified and new ones occur in strong field environ- 
ments [1]. For instance, second-order processes become 
more substantial, e. g., double photon emission [3J [3]. 
Thus, the quantum electrodynamic treatment of such 
processes is necessary when field strength is comparable 
with the critical one {Be = rri^c^/eh « 4.4 • 10^^ G). 

Strong enough constant magnetic field is not feasible in 
laboratory at the present time. Nevertheless, it is possi- 
ble to observe quantum electrodynamic (QED) processes 
in a strong magnetic field in experiments on heavy ion 
collisions. If the impact parameter has order of mag- 
nitude ~ 10~^° sm, then the magnetic field of moving 
ions can approach magnitude of ~ 10""^^ G in the region 
between the ions, while electric fields compensate each 
other. 

At the present time, FAIR (Facility for Antiproton and 
Ion Research) is under construction at the GSI Helmholtz 
Centre for Heavy Ion Research, Darmstadt, Germany. 
One of the goals of the FAIR project is to test QED in 
strong electromagnetic fields. Experiments on observa- 
tion of QED processes in strong magnetic fields in ion 
collisions are possible in the frame of FAIR project. 

Note that the process of pair production by an electron 
has been experimentally observed in an intense laser field 
at SLAC National Accelerator Laboratory 4 . After the 
SLAC experiment, pair creation in laser-proton collisions 
or in countcrpropagating laser beams was studied in a 
number of works, e. g. [5]-|12j. 

Electron-positron pair production by an electron in in- 
tense laser wave was numerically studied in Ref. |13j . 
In particular, the authors considered both resonant and 
nonresonant regimes of the process. In Ref. [T3] the tri- 
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dent pair production amplitude in a strong laser back- 
ground was calculated. 

Pair production by an electron in a magnetic field was 
first studied by T. Erber [15]. In Ref. the rate of 
a cascade of photon emission process followed by pho- 
toproduction in a magnetic field has been estimated for 
both cases of real and virtual intermediate photon. 

In the high-energy limit the considered process in ar- 
bitrary homogeneous constant electromagnetic field has 
been studied in Ref. .16. . 

In Ref. |17) kinematics of the pair production in a mag- 
netic field was considered and the expressions for the to- 
tal process rate containing integrals over orbit centers 
coordinates were obtained. 

The purpose of the present paper is to calculate the in- 
tegrals and obtain the explicit analytical expressions for 
the process rate. The resonant case is studied, when the 
rate factorizes and can be expressed via the product of 
the rates of the corresponding first-order processes. It is 
assumed that all final particles occupy the ground Lan- 
dau level. The explicit analytical expressions for the total 
rate are obtained for subcritical magnetic field strength, 
B<B,. 

Using Nikishov's theorem [THl HH] the obtained result 
has been compared with the experiment on observation 
of pair production by an electron in a laser field 

Relativistic units (fi = c = 1) are used throughout the 
paper. 



II. PROCESS RATE 

Feynman diagrams of the considered process are shown 
in Fig. [ij where the double lines represent the solutions 
of Dirac equations in a magnetic field. 

The process is studied near the threshold, when the fi- 
nal particles occupy the ground Landau level. Lorentz 
transformation does not change magnetic field when 
passing to a reference frame moving along the field. Thus, 
without loss of generality the longitudinal momentum of 
the initial electron can be chosen equal to zero, — 0. 

The corresponding probability amplitude can be writ- 
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Xq, Xoi, and X02 are the x coordinates of the classical or- 
bit centers of the initial and final electrons, respectively, 
and is the x component of intermediate photon mo- 
mentum. 

The purpose of this paper is to carry out integration 
in Eqs. Q - (|6]) and obtain the explicit analytical ex- 
pressions for the process rate in the resonant case. The 
integral D in Eq. ([7| can be expressed in the form 



FIG. 1: Feynman diagrams of the process of e+e -pair pro- 
duction by an electron in a magnetic field. 
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where a is the fine structure constant and Df^^, is the 
photon propagator. 
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where = l\/k\{l — k)\ are binomial coefficients and 
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The integrand has a singularity when the condition 

< is true, and the value of the integral (j?]) is small if 
is positive. Thus, it is necessary to consider the case 

< 0, when the inequality —Vl < s < VL is true. 

The resonant divergence results in the infinite value of 
the process rate. To eliminate the divergence, one should 
introduce a width of the intermediate state A in accor- 
dance with Brcit-Wigner prescription i,20j and replace 



and g^i, is the metric tensor. 

The process rate is defined by the following equation: 



dW 
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Here, S is the normalizing area, d^p = dpydp^. 

The general expressions for the process rate look like 
[Til 



W+ 



37r2^/3^! ' 
W- - bW+ 



(4) 
(5) 



where the superscript denotes initial electron spin projec- 
tion, I is the Landau level number of the initial electron 
and b = B/Bc, B is magnetic field strength and Be is 
critical field strength. The integral Y has the form 



Y = 



ds du 



D 



where 



D 



f {s + iqY ^-^^-2iuQ 
J.2 _ q2 



dq. 



(6) 



(7) 



Here, the following notations are used: 

s = mfl{xQ — Xqi), 

u = mfl{xo - X02), 

q = kx/mV2b, 
r2 = Jl2_^2^ 

^^2 = 2/b. 



(8) 
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A 

mb 



(11) 



The integration in can be carried out analyti- 
cally (see the Appendix) and the result is represented 
by Eq. (lAlOl). 



As noted in the Appendix, the quantity Dq contains a 
divergence in the point s = fl. Thus, when substituting 
Eq. ^ to Eq. ([6|, the summands with fc > 1 can be 
neglected: 



D = 



^ — 2iup 



2ip 



erfc(M — ip) + 

+e^'''PeT{c{-u-ip)} . (12) 



Taking into account, that the width A is small, the 
integration over ds, du in Eq. ^ can be carried out 
analytically too. After the corresponding calculations 



Eqs. (A12)-(A26)] the expression Y takes the form 



Y = bn'^y^ 



n 



2l„-2Q.^ 



T{l + \/2) 



A/m 



(13) 



where T{1 -\- 1/2) is the gamma function. 

Averaging the rate over the initial electron spin pro- 
jection, finally we obtain (in CGS units) 



W 



2 , mc2\ V'i'f^^'e"^^^' + 1/2) 



6^/3 A/m 



(/!)2 



(14) 



The quantity A in Eq. ( |14[ ) should be considered as the 
total width of the intermediate state. The main contri- 
bution to the width is made by the total radiation rate of 
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FIG. 2: The dependence of the total process rate on magnetic 
field strength and the initial Landau level number. 



the initial electron. There are a number of works related 
to this problem, e. g. [2T]-[26]. 



As an example, let us calculate the rate ( 14 1 when field 
strength is 6 = 0.1 {B ^ 4.4 • 10^^ G). In this case the 
threshold Landau level number is / = 40 and 



A«3.9-10i^ (s-i), 
IF =1.2-10'' (s-^). 



(15) 
(16) 



The dependence of the rate (14) on the magnetic field 
strength is shown in Fig. [2j 



III. FACTORIZATION 

One can see that the main contribution to the process 
rate is made by the resonant mode. In this case, the 
total rate factorizes and can be expressed via the prod- 
uct of the rates of the first-order processes of magneto- 
bremsstrahlung and e'^e~ pair production by a single 
photon [211 [27]. With account of the threshold condi- 
tion E « 3r7i, and consequently bl = A and Z 3> 1, the 
following expression can be found: 



W 



3^6 



(17) 



Here, We.^-^e and W^^^e 



+ are the rates of the corre- 
sponding firs-order processes, cyclotron radiation and 
pair photoproduction, respectively: 



>7e 
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be- 
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where 5E = E — 3m and E = m\/l + 2lb is the incident 
electron energy. 

Note that Eq. ( 19 1 does not take account of the state 
widths and diverges if 5E goes to zero. In order that the 
final particles were not allowed to occupy exited energy 
levels the condition 5E < mb should be fulfilled. For 
exapmple, let 5E be ^mb and I = 40, then b = 0.10375 
and Eqs. (18), (19 1 give the following numerical values: 



We-,-te = 2.1 • 10" (s-1), 
VFT,_^ee+ = 7.9 • 10^ (S"l). 



(20) 

(21) 



10^ s-1 while Eq. (14) 



In the review |15j resonant pair production by an ele- 
cron was considered as a cascade of synchrotron emission 
and photoproduction. However, the rates of radiation 
and photoproduction in the high-energy limit were used 
that imply both initial and final states to be ultrarel- 
ativistic. This approach is not applicable near the pro- 
cess threshold when the final particles occupy the ground 
Landau level. 

Moreover, in Ref. [15] the resonant width is not taken 
into account too. In fact, decay time of the virtual state 
assumed to be equal to the half of the time of observation. 

As a result, the approach of Ref. overestimates 
the process rate near the threshold. For the above pa- 
rameters and time of observation equal to the dou bled 
radiative decay time, it yields 4.7- 
gives about 1.2 • lO'^ s^^. 



IV. DISCUSSION 

As stated in the introduction, critical or subcritical 
magnetic field is not feasible in laboratory conditions. 
On the other hand, QED processes have been already 
observed in SLAG experiments involving the interaction 
of an intense laser with an electron beam [H [5S1 [25] . 
In Ref. [4] observation of e+e" pair production by en 
electron in laser field was reported. About 100 positrons 
have been observed in 21 962 collisions of a 46.6 GeV 
electron beam with green (A — 527 nm) terawatt laser 
pulses for which rj = 0.36, where rj — eyJ~AKA^/mc^ and 
Afj_ is four- vector potential of the laser wave. 

The positrons were interpreted as arising from Gomp- 
ton back scattering followed by the multiphoton Breit- 
Wheeler reaction, 



e 



n'ujQ 



+ 6-^ 



(22) 
(23) 



where ujq denotes laser photons. Such a two-step process 
was distinguished from the less probable trident reaction 



n ujQ 



+ e 



(24) 



Nevertheless, it is impossible to observe the interme- 
diate photon Lu' without destroying the whole process. 
The photon should be represented by an internal line in 
the Feynman diagram and by a photon propagator in 
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the probability amplitude (but not by a wave vector). 
Consequently, to develop a consistent theory, one should 



consider the more general trident reaction Eq. ( 24 ) . 

However, when kinematics allows on-shell intermedi- 
ate state (so-called resonance), the Feynman diagram of 



the trident process ( |24[ ) decomposes into two first-order 
diagrams corresponding to the processes (22 1, (23). In 



this case the total rate can be expressed via the rates of 



the processes (22), (23) with some additional coefficient, 



that can be obtained only in the frame of the full theory. 

It is necessary to note that Nikishov and Ritus [TS] 
have proven the form of the expression of the process 
rate to be the same for any external field if the rate is 
expressed in terms of gauge invariants and velocity of the 
incident particle is ultrarelativistic. In Ref. [IH] the rates 
of one-vertex processes were obtained in the case of a 
laser field. If the variability of the laser field is irrelevant, 
the obtained expressions reduce to the rate of the pro- 
cesses in crossed electric and magnetic fields, when £ _L _B 
and 8 — B. The total rates of such processes depends on 
the single invariant parameter e'^(F^^p^Y' /m^ where F^jy 
is the electromagnetic tensor and pi, is 4-momentum. It 
allows to pass to the general case of arbitrary constant 
field. In this case rates depend also on two other param- 
eters e^F'^^/m'^ and ie^e^^yXcrF^^^ F'^'^ (they are equal to 

zero \i £ ^- B and £ = B). 

However, since feasible fields are much less than the 
critical one w? /e, these additional parameters are much 
less than unity. On the other hand, if the particle energy 
is high enough, then these parameters are much less than 
the first one as well and could be omitted. Therefore, the 
obtained rates are applicable in the case of arbitrary con- 
stant field, if the incident particle has relativistic energy. 

In particular, considering F^^^^ as a magnetic field, Nik- 
ishov and Ritus have obtained the results of Klepikov [55] 
for intensity of a photon emission by an electron and for 
the rate of pair production by a photon in a magnetic 
field. 

The physical reason is that due to Lorentz transforma- 
tion arbitrary electromagnetic field goes to almost equal 
and almost perpendicular electric and magnetic fields 
when passing to the rest frame of the relativistic particle. 

Thus, it is possible to compare the analytical result for 
the case of magnetic field with the experimental data of 
Ref. [4]. 

If a relativistic electron propagates opposite to elec- 
tromagnetic wave of field strength fi, then it experi- 
ences the field strength of £q = 2^8 l in the rest frame, 
where 7 is the gamma factor. On the other hand, if 
an electron moves perpendicular to a magnetic field Beg, 
then the field strength in the rest frame is approximately 
£oeq = iBeq- Comparing £q and foeq one can see that 
strength of the equivalent magnetic field in the lab frame 
is 



Beq — 2£i 



(25) 



field should take into account both electric and magnetic 
fields of the electromagnetic wave. 

In order to pass to the case of alternating field of an 
electromagnetic wave, the rate (14) for the process in 



a magnetic field should be averaged over the wave period 
to obtain the equivalent process rate in laser field W^q 



Weq^^ J WiBeqSin. 



(26) 



Equation ( 26 ) allows us to compare the rates of processes 



in a magnetic field and in an intense laser wave. 

However, Eq. ( 14 ) is true near the process threshold 
only, when the condition E « 3m is fulfilled. There- 
fore, it is necessary to calculate the rate in the moving 
"threshold" frame where the electron energy is equal to 
E w 3™, and threshold conditions are fulfilled explic- 
itly. The amplitude value of equal magnetic field in the 
threshold frame is Beq w 6.1 • 10^^ G and, consequently, 
6 « 0.14. 

It should be noted that in the SLAC experiment pair 
production has been observed near the threshold too ||] . 
Although the electron beam energy was 46.6 GeV, the 
major part of this energy was the energy of rectilinear 
motion of the mass center. 

It is possible to estimate the electron-laser interaction 
time in the laboratory frame AIl and the number of 
electrons in the interaction region Ni„t using the data 
from Ref. [3]: the electron beam size is ^ 25 x 40 /xm^, 
bunches contained ~ 7 • 10^ electrons, laser beam focal 
area is 30 ^m^, beams crossing angle is 17°. Thus, AIl ~ 
50 fs, N,„t - 2.8 • 10^ 



Note that to calculate the rate Eq. ( 14 ) it is necessary 



to take into account limited interaction time as well as the 
radiative width (15). Therefore, the intermediate state 



width is a sum of the radiative width and the quantity 
1/AtT where Atr — Ati^j^ is laser-electron interaction 
time in the threshold frame. 

The number of produced pairs can be estimated ac- 
cording to the expression 



(27) 



Note that factor 2 arises because equivalent magnetic 



where fc = 21 962 is the number of collisions of the elec- 
tron and laser beams [3|. 

The corresponding value of ~ 80 events is in reason- 
able agreement with the experimental result of 106 ± 14 
indicated in Ref. \p . 

Note that the authors of Ref. |1] pointed out the possi- 
ble residual background of about 2 x lO^'' positrons/laser 
shot due to interactions of Compton backscattered pho- 
tons with beam gas. If the data are restricted to events 
with 77 > 0.216, one can find 69 ± 9 positrons, and the 
agreement of their number with theoretical estimations 
is improved. 

Thus, in the present work the analytical expression 
for the rate of electron-positron pair production by an 
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electron in a magnetic field near the process threshold 
was obtained. The number of e^e~ pairs created in the 
SLAC experiment was estimated using Nikishov's theo- 
rem. The obtained value is in reasonable agreement with 
experimental results as well as with numerical calculation 
oftheRef. [13|. 

We thank V. Yu. Storizhko and S. P. Roshchupkin for 
useful discussions. 



Appendix A: Calculation the integrals 

To take an integral of the form 



-dq 



(Al) 



it is convenient to use the apparent relation 

— — q^ 





(A2) 



The integral ( Al I takes on the form 



= e 



p^ — (f- 



dq+e 



""e'p'-^dt. (A3) 



The first integral in (A3) can be found using Jordan's 
lemma. The result is 



2iuq 

dq = Le^^Hi, 

ip 



p^ — cp 



(A4) 



To find the second integral one should use the substitu- 
tions 



(T+ = p^ft + i\u\l\Jt, 
(J- = p\/t — i\u\l 



(A5) 



After simple calculations the result of integration takes 
on the form 



TT 

2ip 



3-2vl«lerfc(|u| - ip) - e^^^l^lerfcditl -f ip)] . (A6) 



Finally, substituting Eqs. ( A4 1, ( A6 1 to Eq. ( A3 ) the re- 



sult for Z?o can be expressed as 

Do = — re-2v«erfc(M - ip) + e^^^erfcf-w - ip)] . 

zip 

(A7) 

The above expression is valid for both u > and u < 
cases. 

Integrals containing q'' can be reduced to the consid- 
ered one using the derivative with respect to the param- 
eter u: 



Taking into account the relation 

Hn{x) = (-l)"e" -^e- 



(A9) 



where Hn{x) is Hermite polynomial, the explicit form of 
Dk can be expressed as 



Dk = ^—p'' [e-2»Perfc(^x-zp)- 

Zip 

+ (-l)'=e2*"''erfc(-u- 



+ E C^ni^^P)'-"" iH,n-l{u ip) + 

+ . (AlO) 

Note that the value Dq is in inverse proportion to p 
and contains a divergence in the point s = fl. On the 
contrary, the valu e Dh is finite for k ^ 1. Indeed, the 
summands in Eq. (AlO) contain factors p'^^^ and p™-*^-! 



and apparently do not diverge when k > m. When the 
conditions m — k and p = arc true, then the second 
summand contains the expression 



[Hk-i{u) + [-l)''Hk-i{-u)] =0 



(All) 



where the relation Hn{—x) = {—l)'^Hn{x) is used. 

Let us proceed to calculating the integral Y in E g. d6| . 
Taking into account, that the quantity D in Eq. Jl2| is 
an even function of b, the integral over db in Eq. ([6|) can 
be expressed as 

\e-' Dfdu^ .\. [J1 + J2] (A12) 



where 



Ji - 



|e~2»Perfc(u-ip)pdu, (A13) 



''erfc(u — ir)crfc(— u + ir)du. (A14) 



After integration by parts the quantity Ji takes the 
form 



Ji 



1 



3{P) 



y/n'^p 
e-(«+v)%rfc(w- ip)dw. 



(A15) 
(A16) 



The parameter p can be eliminated from the argument of 
the exponent by introducing the new variable t = u + ip: 



Jip) 



erfc(t — 2ip)dt. 



(A17) 



The derivative of j{p) with respect to p reduces to the 
Poisson integral and takes the form 



j'{p) = 2V2ie^p\ 



(A18) 
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This differential equation can be easily solved. Finally, Taking into account that 



after substituting the result into Eq. (A15) the quantity 
Ji takes on the form 

jR[e-^'gerfc(-zp^)] 

W) ■ ^ ^ 

Note that Ji can be expressed as 

Jr « ^ + (A20) 



The integral J2 in Eq. (A14| can be calculated in the 
same way and looks like 



erf(^^r^/2) 



(A21) 



^<P) = ^47^+? - (A23) 



and introducing a new variable x — s/Q, the quantity Y 
can be transformed to 
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where 6 — .g^/il. When S goes to zero the integral over 
X in the above expression converges to 



The value of the integral over s is determined by the 
region in the vicinity of the point s = due to the p res- 
ence of the factor s^' in the integrand in Eq. (A12). In 
the points s — zLil the first summand in Ji Eq. (A20) 



goes to y 2/5, while the second one and the quantity J2 
go to ±^/8/^. Thus, 



Y = TT^e 



n ^21 



ds 



(A22) 



Thus, 



1 r(i/2)r(f + i/2) 
V2 ni + i) 



-2n- 1 ni + 1/2) 
g II ■ 



(A25) 



(A26) 
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